Chern Numbers for Spin Models of Transition Metal Nanomagnets 
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We argue that ferromagnetic transition metal nanoparticles with fewer than approximately 100 
atoms can be described by an effective Hamiltonian with a single giant spin degree of freedom. The 
total spin S of the effective Hamiltonian is specified by a Berry curvature Chern number that char- 
acterizes the topologically non-trivial dependence of a nanoparticle's many-electron wavefunction on 
magnetization orientation. The Berry curvatures and associated Chern numbers have a complex de- 
pendence on spin-orbit coupling in the nanoparticle and influence the semiclassical Landau-Liftshitz 
equations that describe magnetization orientation dynamics. 
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Both molecular nanomagnetshj, JjA |3j and ferromag- 
netic transition metal clustersjj, 0, L| have been ac- 
tively studied over the past decade. Interest in molec- 
ular nanomagnets has been due primarily to their po- 
sition near the borderline between quantum and clas- 
sical behaviors 0, 0, @- For ferromagnetic transition 
metal nanoparticles, on the other hand, interest has been 
spurred mainly by classical physics issues relevant to in- 
formation storage [ij [ToL ITl| and by the interplay between 
collective and quasiparticle degrees of freedom[r3. Ilflj . 
The present work is motivated by the observation that 
the low energy physics of small transition metal clusters 
can be described by an effective Hamiltonian with a sin- 
gle giant spin degree of freedom, like that of a molecu- 
lar magnet. A transition metal nanoparticle will behave 
like a molecular magnet when the energy scale associ- 
ated with its collective magnetization orientation, the 
anisotropy energy, does not exceed the smallest energy 
scale associated with its quasiparticle degrees-of-freedom, 
the single-particle level spacing 5. When bulk density-of- 
states and anisotropy energy values are used to estimate 
the particle-size at which this condition is satisfied, the 
cubic transition metal ferromagnets Fe and Ni are pre- 
dicted to act like molecular magnets when the number of 
atoms Na is smaller than ~ 1000, while Co is predicted 
to act like a molecular magnet for Na smaller than ~ 100. 
(See Table [J. It is remarkable that an extensive (oc Na) 
energy scale, like the total anisotropy energy, is smaller 
than a microscopic energy scale (oc N^ 1 ), like the level 
spacing, at a relatively large particle number. This sur- 
prising property arises from a combination of relatively 
weak spin-orbit coupling in the 3d transition-metal se- 
ries and the itinerant character of transition metal ferro- 
magnetism. The magnetocrystalline anisotropy energy in 
these systems in the bulk is five to six orders of magnitude 
smaller than the magnetic condensation energy. Even 
accounting for the substantially larger anisotropy energy 
per volume expected in typically shaped nanoparticles, it 
seems clear that the total width of the anisotropy energy 
landscape will often be smaller than the level spacing for 
Na smaller than ~ 100. The physics that controls the 
Hamiltonian and, most centrally, the total spin of these 
malleable molecular magnets is the subject of this Letter. 



TABLE I: 8 Na is the bulk mean-level spacing times the 
number of atoms. Eb is the bulk coherent rotation anisotropy 
energy barrier (For cubic systems Eb/vo\ is one third of mag- 
netic anisotropy constant Ki, whereas in the uniaxial case it 
is equal to the magnetic anisotropy constant.) NX is the num- 
ber of atoms at which 5 and Eb become equal. 
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Fe 


695 


0.016 


.00118 


767 


Co 


581 


0.53 


.0367 


126 


Ni 


495 


0.002 


.000114 


2080 



"From Ref. [3 

''Bulk value at room temperature from R.ef, lift 



In the absence of spin-orbit and dipole-dipole interac- 
tions, the collective physics of a transition metal nanopar- 
ticle is easily understood [L4J. Quasiparticle states, cal- 
culated for a particular magnetization orientation (us- 
ing the spin-density-functional theory recipe for exam- 
ple) have definite spin-orientations. The many-particle 
ground state has a total spin quantum number S — 
(N a — N )/2, where N a and N a are respectively the 
number of quasiparticle states whose individual spins are 
aligned with and opposed to the total spin. The 2S + 1 
fold ground-state degeneracy is lifted by an arbitrarily 
weak external magnetic field, selecting a ground state 
with magnetic moment 2SfiB- This simplicity dissolves 
into a complex muddle when spin-orbit interactions are 
included. Individual quasiparticle states no-longer have 
definite spin orientations. It is not immediately obvious 
how to derive an effective Hamiltonian which describes 
the low-energy collective physics, or even how to deter- 
mine the dimension of the quantum Hilbert space for 
the coherent magnetization degree of freedom. We ad- 
dress this issue here using a path integral approach. We 
find that the total spin S of the molecular magnet is 
specified by a Chern number defined by the dependence 
of the many-particle wavefunction phase on magnetiza- 
tion orientation, and propose a procedure for extracting a 
quantum Hamiltonian from the dependence of energy and 
Berry curvature on orientation. Our analysis is based on 
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an approximate imaginary-time quantum action with a 
single magnetization-orientation degree of freedom, n(r): 



Dh[n\ 



EE / dT 



Vn^[n}) ■ — +E[h] 



(1) 



For simple model Hamiltonians in which the interaction 
terms can be represented exactly by auxiliary field func- 
tion integrals 17|, actions of this form can be derived fully 
microscopically 18J. For the realistic description of real 
nanoparticles, however, this action is approximate and 
should be constructed by using the family of spin-density- 
functional theory constrained Kohn-Sham states with net 
magnetization orientation h for an approximate iden- 
tity resolution at energy scales below the single-particle 
level spacing. In the absence of spin-orbit interactions, 
spin-density-functional theory describes the ground state 
in terms of Kohn-Sham quasiparticles that experience 
a scalar potential and an exchange-correlation effective 
magnetic field 19], both of which have a complex de- 
pendence on spatial coordinate. The oricntational de- 
gree of freedom we have in mind for the action is the 
direction n in spin-space of the spin-de nsity - functional 
theory exchange-correlation effective field 20] , or its spa- 
tial average if low-energy states have non-collinear mag- 
netization. In Eq.[^ I 1 ! 7 !"-]) is the Kohn-Sham single- 
Slater determinant state defined by this orientation, 
with spin-orbit terms included in the Hamiltonian and 
E[n] is the density-functional-theory energy. The first 
term in this action is a Berry phase contribution|2l) 
whose role in quantizing small-amplitude magnetic orien- 
tation fluctuations (spin-waves) of bulk transition metal 
ferromagnets has been successfully exploited by Niu 
et. al. |22l 123. |24| . and other authors |25|. When spin-orbit 
interactions are included, both Berry phase and energy 
function E[n] terms will have a non-trivial dependence 
on orientation^,!^. We expect that quantitatively ac- 
curate actions can be constructed in this way for com- 
pletely specified nanoparticles; in the rest of this paper we 
discuss some qualitative properties using a tight-binding 
model^. 

It is useful to express the Berry phase in a gauge invari- 
ant form|2ll l28j . During its imaginary-time evolution, 
the unit vector n traces a closed trajectory on the unit 
sphere. By Stokes' theorem we can rewrite the closed 
line integral and express the Berry phase in terms of a 
surface integral over the area enclosed by the path 

S Bc rryN = f dn-/*|V ft *\ = / V ft x faf\ V fi *Yn da 

J J area 

(2) 

The Berry curvature 



C[h] =iV fl x (*[n] S7nV[n 



(3) 



is gauge invariant and the Chern number S, defined as 
the average of the curvature over the unit sphere, is re- 
quired to be half of an integer |2j|. The Chern number 
appears below as the total spin in the effective model for 
the nanoparticle, hence the notation chosen above. 



To discuss the Chern numbers, it is useful to start 
from the case of no spin-orbit coupling, for which S = 
(N a — N )/2. To see this, consider the Kohn-Sham single- 
Slater determinant many-particle state which describes 
magnetization in direction n: 

= J] [ U (n)cj T +u(n)c!j J] [-«*(n)c} T +u(n)ctJ |0) , 



where cl 



(4) 

s are electron creation operators defined 
with the arbitrary polar angle of the coordinate sys- 
tem used to measure magnetization orientations as the 
spin-quantization axis, and i and j are orbital indices 
running respectively over N a majority-spin orbitals and 
N Q minority-spin orbitals. We have made the gauge 
choice u(h) = cos(0/2), v[n) = e~ lv sm(9/2) for the 
S = 1/2 quasiparticle coherent states, where 9 and 
ip are the polar and azimuthal angles used to specify 
n(9,ip). Given this state, the Berry connection |2ll l2q. 

A[n] ee i/^[n] Vft^fn] ), is easily calculated 



. r du ^ 4 dv ' ^ . r dv* ^ du 
J L dn dni ^— J L dn dn 

i 



1 - cos(0) „ _ 



dni V dn dn 

j 

N a -N l- cos(9) 



sm(9) 



0,(5) 



It follows that the Berry curvature C[n] = x A[h] is 
constant and equal to 



C = 



(6) 



In this case each majority-spin Kohn-Sham quasipar- 
ticle contributes 1/2 to the Chern number while each 
minority-spin quasiparticle contributes —1/2. 

When the spin-orbit interaction is present, the Berry 
curvature is no longer a constant, but the Chern number 
S must still be half-integer. The total Chern number S of 
the nanoparticle can still be considered as the sum over 
contributions from each quasiparticle state .22]. Individ- 
ual quasiparticle Chern number properties are familiar 
from related studies, especially in connection with the 
topological interpretation of the quantized Hall conduc- 
tance in a two-dimensional electron gas[2j| 13(3, OH], 132] . 
They can change only when quasiparticle level degen- 
eracies occur. In nanoparticles without special symme- 
tries, level degeneracies will occur at discrete values of 
h, and at isolated values of any microscopic Hamiltonian 
parameter. When two levels intersect and then come 
apart again as a parameter is varied, then their individ- 
ual Chern numbers may not be conserved but their sum 
is conserved |3l|. As we now show, spin-orbit interac- 
tions lead to surprisingly ubiquitous and large deviations 
of quasiparticle Chern number values from their ±1/2 
values in the absence of spin-orbit coupling. 

In Fig. n we plot the variation of quasiparticle Chern 
numbers for representative groups of contiguous levels of 
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FIG. 1: Chern numbers of quasiparticle energy levels of a 
25-atom Cobalt nanoparticle, as a function of the spin-orbit 
coupling strength £. The nanoparticle is modeled by a tight- 
binding Hamiltonian|27j| . In each panel I, II, (and sometimes 
III) label two (or three) contiguous orbitals whose energies 
are ordered in ascending order. Deviations of individual 5* 
from their zero spin-orbit coupling values ±1/2 occur when 
two levels cross. Note that the sum of the two Chern numbers 
is conserved at crossing. 
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FIG. 2: Chern numbers of all the individual quasiparticle 
energy levels of a 25-atom nanoparticle for a given value of 
£. As a result of repeated level crossings, several of these 
numbers differ strongly from ±1/2. 



a nanoparticle modeled by a tight-binding model |27| , as 
a function of the spin-orbit coupling strength £. Typi- 
cally we find that individual Chern numbers of pairs of 
levels undergoing level crossing experience a change of 
±1, with their sum remaining unchanged. These changes 
can be thought of as representing changes in the orbital 




FIG. 3: Planar projection of the total (negative) Berry cur- 
vature —C[n] for a 25-atom nanoparticle. The strength of the 
spin-orbit coupling £ = 165 meV is chosen in proximity to a 
value at which the Fermi level and the first unoccupied level 
cross. Level crossing occurs for a given h and its opposite —n. 
In correspondence of these two values, — C[h] has large peaks. 



contribution to the effective angular momentum of indi- 
vidual orbitals. Note that level crossings always occur 
in pairs, for a given value of n and its opposite — n|33|. 
As a result of repeated level crossing, some of the indi- 
vidual Chern numbers end up with values very different 
from ±1/2 as shown in Fig. [3 where we plot all individ- 
ual quasiparticle Chern numbers for a given £. If a level 
crossing occurs between the Fermi level and the first un- 
occupied level, the total Chern number will change from 
the original zero spin-orbit value (N a — N a )/2. In the 
vicinity of such a level crossing, the Berry curvature C[fi] 
will deviate strongly from its average value as a function 
of [h]. An example of this is shown in Fig. |3 where we 
plot C[n] at a value of £ for which the Fermi level and 
the first unoccupied level are almost degenerate in two 
directions of h and the Berry curvature has sharp peaks. 

Assuming only that the Berry curvature is positive- 
definite, our low-energy effective action can be mapped 
to that of a quantum spin Hamiltonian by making a 
change of variables that transforms the Berry curvature 
to a constant. The Hamiltonian representation is easier 
to use for explicit calculations of collective tunneling am- 
plitudes, non- linear response to external electromagnetic 
fields and other relevant properties. To be explicit, we 
change variables from u = cos{9) and tp to u' and ip' 
defined by 

u = -1+2^5 J du " J dp"C{u",p") (8) 
With this change of variables C(u, tp) dudip = S du'dip' 
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and the real-time action for a path can be written 



dt ' (^'^mkm^^ 



(9) 

where A = S(l - cos(0')) / sin(6')4>'. This is the quan- 
tum action for a total spin quantum number S with the 
classical Hamiltonian\2£fy specified by the second term of 
the integrand, i.e. 



H[n'(t')] = {S,h'(t')\H\S,n'(t')) 



SE[n(h'(t'))] 
C[n(h'(t'))] 



where TL is the quantum Hamiltonian of the spin system 
and \S, n'(t')) is a spin-S* coherent state parametrized by 
the unit vector n'(t'). Given the energy and Berry cur- 
vature functions, this quantum Hamiltonian can always 
be explicitly constructed. Applications of this procedure 
will be presented elsewhere. 

A non-constant Berry curvature affects the Landau- 
Lifshitz equations describing the precession motion of 
the nanoparticle magnetic moments. These equations 
are equivalent to the Euler-Lagrange equations of mo- 
tion derived from the real-time action <S sp i n [n'] in the 
semiclassical (large S) approximation: 



h cl (t') = n c i(t') x 



d[H[h'(t')]/S 



dh'(t') 



(11) 



where n c \(t') is determined by the saddle-point equations 

SS S y in [n'(t')] 



&h'(f) 



«cl 



0. Using the expression of H in Eg. 1111 

we can see that a non-constant Berry curvature modifies 
the precession rate of the magnetic moment fluctuations. 
This effect is particularly important when the Fermi level 
and the first unoccupied state are close to a degeneracy 
point for some value of n, since as shown in Fig El C\p\ 
can then deviate considerably from the Chern number S. 



In conclusion, we have derived an effective spin Hamil- 
tonian with a single giant spin degree of freedom that de- 
scribes the low-energy physics of a small metallic nano- 
magnet. The dimension of the Hilbert space of the ef- 
fective Hamiltonian is given by a Berry curvature Chern 
number which has a non-trivial dependence on spin-orbit 
coupling strength and nanoparticle details. We would 
like to thank G. Canright and W. Wersdorfer for help- 
ful conversations. This work was supported in part by 
the Swedish Research Council under Grant No:621-2001- 
2357, by the faculty of natural sciences of Kalmar Uni- 
versity, and by the National Science Foundation under 
Grants DMR 0115947 and DMR 0210383. 
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